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IMAGES OF AN EQUATORIAL OUTFLOW IN SS433 
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ABSTRACT 

We have imaged the X-ray binary SS 433 with unprecedented Fourier-plane coverage at 6 cm using 
simultaneously the VLBA, MERLIN, and the VLA, and also at 20 cm with the VLBA. At both wave- 
lengths we have securely detected smooth, low-surface brightness emission having the appearance of a 
'ruff' or collar attached perpendicularly to the well-studied knotty jets in this system, extending over at 
least a few hundred AU. We interpret this smooth emission as a wind-like outflow from the binary, and 
discuss its implications for the present evolutionary stage of this system. 

Subject headings: Stars: Binaries: Close, Radio Continuum: Stars, Stars: Individual: SS433 



1. INTRODUCTION 

SS 433 is famous as the first known relativistic jet source 
in our Galaxy. Red- and blue-shifted optical lines, indicat- 
ing velocities up to 0.26c, were discovered by Margon et al. 
(1979a, b) and interpreted as gas accelerated by oppositely- 
directed jets (Fabian & Rees 1979; Milgrom 1979). Margon 
(1984) successfully fit a kinematic precessing jet model, 
finding an instrinsic jet speed of rj 0.26c, a precession pe- 
riod of 163 days, a cone opening angle of ~ 20°, and an 
inclination to the line-of-sight of ~ 80°. This model was 
spectacularly confirmed by subsequent radio imaging (e.g., 
Hjellming & Johnston 1981; Vermeulen et al. 1987; Fejes, 
Schilizzi, & Vermeulen 1988; Vermeulen et al. 1993), which 
showed precessing twin jets with structure on scales from 
milli-arcseconds to arcseconds. X-ray emission lines were 
also discovered which mimic the behavior of the 'moving' 
optical lines (Kotani et al. 1994; Marshall et al. 2001). On 
the basis of these and other photometric data, SS 433 is 
believed to be a binary consisting of a compact object [a 
black hole (Zwitter & Calvani 1989; Fabrika & Bychkova 
1990) or a neutron star (Filippenko et al. 1988; D'Odorico 
et al. 1991)] and an O- or B-type star (Margon 1984). Dis- 
tance estimates range from 3.1 kpc (Dubner et al. 1998) 
to 5.5 kpc (Hjellming & Johnston 1981). Here we report 
the first results from simultaneous radio observations with 
the Very Long Baseline Array (VLBA), the Multi-Element 
Radio Linked Interferometer Network (MERLIN), and the 
Very Large Array (VLA). In addition to the well-known 
relativistic jets, which will be examined in a series of com- 
panion papers, these data reveal smooth radio emission 
centered on the radio core, extending in a huge 'ruff' or 
collar up to > 40 milli-arcsec either side of the center. 

2. OBSERVATIONS AND IMAGES 

We observed SS433 on 6, 7, and 8 March 1998 for four 
hours each day with 26 antennas of the VLA; for lOhrs 
each day with the 6- antenna MERLIN array; and for 12 hrs 
each day with the 10-antenna VLBA together with one 
VLA antenna. The VLA interleaved observations at ~ 20, 
6, and 2 cm; the VLBA, at ~ 20, 13, and 6 cm; and MER- 
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LIN observed only at ~ 6 cm. The images shown here 
were made from the 7 March 1998 data, combining all the 
arrays at ~ 6 cm, and using the VLBA alone at ~ 20 cm. 
Images from the independent data taken on 6 and 8 March 
1998 confirm all the features reported here. 

The 6 cm image shown in Figure 1 was made by com- 
bining the 7 March 1998 data from all three arrays. MER- 
LIN observed at 4.993 GHz with a bandwidth of 15 MHz 
in each of the two independent circular polarizations, us- 
ing J 1907+0127 as a phase-reference source. The VLBA 
observed at four independent frequency bands (IF pairs), 
each 8 MHz wide, in each of the two circular polarizations. 
Two IF pairs (centered on 4.990 and 4.998 GHz) directly 
overlap MERLIN's frequency coverage; the other two are 
not considered here. J 1907+0127 and J 1929+0507 were 
used as phase-reference sources, sandwiching 120 s scans 
on SS 433 between 70 s scans on those two calibrators. The 
total time on SS 433 (after flagging) at this wavelength was 
2.7 hrs. The VLA observed with two IF pairs, each 50 MHz 
wide, in each of the two circular polarizations, centered 
on 4985.1 and 4614.9MHz. The calibrator 1950+081 was 
used for phase referencing, while occasional observations 
of 3C 286 set the fundamental flux scale. Initial phase and 
flux calibration, and fringe-finding for the VLBA, were 
carried out using standard procedures in NRAO's Astro- 
nomical Imaging Processing System software. The (in- 
dependently derived) flux scales of the three instruments 
agreed to better than a few per cent, as determined in over- 
lapping sections of the uv-plane. Finally, the MERLIN 
data, the overlapping IF pairs from the VLBA, and the 
higher-frequency IF pair from the VLA were concatenated 
to produce a final, combined data set with full sensitivity 
on spatial scales ranging from ~ 10~ 3 to ~ 10 arcseconds. 

The VLBA observed SS433 using four IF pairs, each 
covering 8 MHz in each of the two independent circular 
polarizations. The central frequencies were 1.404, 1.505, 
1.651, and 1.658 GHz, chosen to be spread as far apart as 
possible within the 18-21 cm band. The same observing 
scheme was used as at 6 cm. The image shown here how- 
ever is referenced to J 1928+0507. To avoid uncertainties 
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associated with a spatially- variable spectral index, each IF 
pair was independently imaged (e.g. Fig. 2). 




19 11 49.572 49.570 49.568 49.566 49.564 49.562 49.560 
RIGHT ASCENSION (J2000) 
Grey scale flux ranges -1 .26 93.42 MilliJY BEAM 
Cont peak flux = 9.341 7E-02 JY/BEAM 

Fig. 1.— 4.99 GHz image of SS433 on 7 March 1998, com- 
bining VLBA, MERLIN, and VLA data. The restoring beam is 
9.5 X 8.2 milliarcscconds (FWHM), oriented at a position angle of 
— 18?1. Contours are 0.5 X \/2 n mjy/beam, n= 1,2, ... 

10 20 30 40 50 




19 11 49.576 49.574 49.572 49.570 49.568 49.566 49.564 49.562 49.560 49.558 49.556 
RIGHT ASCENSION (J2000) 
Grey scale flux range= 0.50 54.00 MilliJY BEAM 
Cont peak flux = 4.6058E-02 JY/BEAM 

Fig. 2.— 1.6 GHz VLBA image of SS433 on 7 March 1998. The 
restoring beam is 13.8 X 7.5 milliarcscconds (FWHM), oriented at a 
position angle of — 6?7. Contours are 0.5 X \/2™ mjy/beam, n= 1, 
2, 13. 



3. A THERMAL 'RUFF' ABOUT THE CENTER OF SS 433 

Apart from the well-known knotty jet, the most striking 
feature of these new images is the smooth emission seen 
perpendicular to the jet at its center. Since this emission 
appears to encircle the main jet, we refer to it as a 'ruff'. 

3.1. Is the 'Ruff' Real? 

Careful testing persuades us that the 'ruff emission is 
real, and not some calibration or imaging artifact: self- 
calibration (avoiding the ruff emission itself) and careful 
data editing cannot remove the emission, which is clearly 
visible after a single round of phase-only self-calibration 
based on the brightest portions of the jet. The emission is 
seen at consistent flux levels, positions, and orientations: 
in MERLIN and VLBA data reduced independently, and 
in data from three consecutive days reduced separately. 



The 'ruff' appears with similar size, brightness asymme- 
try, position, and flux density (sec below), in images made 
in the individual IF pairs covering 18-21 cm and 6 cm. 

3.2. Previous Observations 

Paragi and co-workers were the first to report firm ev- 
idence for equatorial emission distinct from the jet itself 
(Paragi et al. 1998, 1999a,b), from VLBA data taken in 
May 1995. There are some significant differences in what 
we see, for example they detect equatorial emission which 
appears to be rather more blobby than smooth; and while 
our images suggest a 'halo' which connects smoothly to 
(or around) the jet, their naturally-weighted 1.6 GHz im- 
age (fig. 1 in Paragi et al. (1999a)) shows a gap of more 
than 30 mas between the two. It may be that these differ- 
ences indicate genuine variability in this structure. 

Our discovery is the first reported detection of such a 
smooth feature in SS433 extended perpendicular to the 
jets. In part this is because of the high-quality Fourier 
sampling of these observations, which gives good sensitiv- 
ity over a wide-range of spatial scales, and also provides 
much higher fidelity images, less vulnerable to calibration 
and imaging artefacts. 

Evidence of extended emission on much larger angular 
scales may be seen approximately north and south of the 
central circular blob in the 330 MHz image of Dubner et 
al. (1998, figs la & lb). 

3.3. Characterizing the 'Ruff' Emission 

In order to determine the spectral shape of the 'ruff' 
emission we chose wtJ-weighting schemes which gave 
roughly the same dirty beams for both the 6 and the 
18 cm data, CLEANed the images to somewhat below 
the rms noise levels, and finally convolved the resulting 
(model+residual) images to a common 10 x 10 mas Gaus- 
sian beam. The resulting total flux densities, measured in 
identical boxes in all images, which were chosen to avoid 
the jet but include the full 'ruff' emission, are shown in Fig- 
ure 3a. Note that these are lower limits to the total 'ruff' 
emission, since we are excluding any such emission which 
overlaps with the jet at this spatial resolution. The spec- 
tral index for the combined (northcrn+southern) emission 
is a = —0.12 ±0.02 (S v oc v a , where S v is the flux density 
at frequency v) . This is an interesting result since most re- 
solved synchrotron sources are characterized by a. < —0.4; 
indeed, a = —0.1 is normally considered the signature 
of thermal bremsstrahlung emission as is often observed 
in outflows from symbiotic binaries (Scaquist, Taylor, & 
Button 1984; Mikolajewska & Ivison 2001). The compli- 
cation here is that the peak surface brightness corresponds 
to a brightness temperature of (2 — 4) x 10 7 K at 18 cm, 
implying a similar lower limit to the physical temperature 
of a thermally-emitting plasma. 

The distribution of the flux density perpendicular to the 
jet is shown in Figure 36, which suggests that the spectral 
index is indeed almost flat throughout the ruff, and fur- 
ther shows that the emission extends to > 40 mas at our 
sensitivity, or ~ 120 (d/3kpc) AU . Note also that the ruff 
is roughly symmetric about the jet. 

3.4. Implications for X-ray emission 

If the radio emission is bremsstrahlung from a thermal 
population of particles then that same population should 
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co-spatially emit X-rays. The X-ray luminosity is given 
by L x = L ra d x enp(—hv K /kT) where T is the tempera- 
ture of the particles (approximated to be the radio bright- 
ness temperature) , k and h are the Boltzmann and Planck 
constants, v x is the lower frequency of the X-rays whose 
luminosity L x is predicted from the radio luminosity L ra d- 
We measure a brightness temperature Tb>10 7 K which 
over-predicts the X-ray luminosity compared with that ob- 
served (Margon 1984). Allowing for a significant fraction 
of the observed X-rays being emitted by the jets means the 
discrepancy could be as much as an order of magnitude. 
This could in principle be due to the presence of neutral 
material in the vicinity of the ruff which would absorb X- 
rays but not radio emission. Alternatively, it may suggest 
that the particle population is not Maxwcllian. 

4. THE ORIGIN OF THE SMOOTH EMISSION 

4.1. Theoretical background 

The most straightforward interpretation of the radio 
emission is that it arises from mass outflow from the bi- 
nary system that is enhanced towards the orbital plane. 
Such mass loss could cither (i) come from the companion 
(most likely an O or B star), (ii) be a disk wind from the 
outer parts of the accretion disk or (iii) arise from mass 
loss from a proto-common envelope surrounding the bi- 
nary components. The detection of this mass loss may 
have rather important implications for our understanding 
of the evolutionary state of this unique system. It has been 
a long-standing puzzle how SS 433 can survive so long in 
a phase of extreme mass transfer (M > 10 _5 M Q yr _1 ) 
without entering into a common envelope phase where 
the compact object spirals completely into the massive 
companion (for a recent discussion see King, Taam & 
Bcgclman 2000). Since the theoretically predicted mass- 
transfer rate exceeds even the estimated mass-loss rate 
in the jets (M ~ lO^Mgyr" 1 ; Begelman et al. 1980), 
King et al. (2000) proposed that most of this transferred 
mass is lost from the system in a radiation-pressure driven 
wind from the outer parts of the accretion disk (see also 
King & Begelman 1999). A related problem exists in 
some intermediate-mass X-ray binaries (IMXBs). Mod- 
els of the IMXB Cyg X-2 (King & Ritter 1999; Podsiad- 
lowski & Rappaport 2000; Kolb et al. 2000; Tauris, van 
den Heuvel, & Savonije 2000) show that the system must 
have passed through a phase where the mass-transfer rate 
was ~ 10~ 5 Mq yr -1 , exceeding the Eddington luminosity 
of the accreting star by many orders of magnitude, without 
entering into a common-envelope phase, and where almost 
all the transferred mass must have been lost from the sys- 
tem. The observed emission in SS 433 presented here may 
provide direct evidence of how such mass loss takes place. 

4.2. Evidence for mass out-flows 

The existence of a disk-like outflow was first postulated 
by Zwitter, Calvani, & D'Odorico (1991) to explain the 
variation with precession phase of the secondary minimum 
in the photometric light curve. Fabrika (1993) proposed 
a disk-like expanding envelope caused by mass-loss from 
the outer Lagrangian point L2 to explain the blue-shifted 
absorption lines of HI, He I and Fell (see also Mammano, 
Ciatti & Vittone (1980), whose spectrum shows that all 



the emission lines seen in SS 433 have P-Cygni profiles in- 
dicating the presence of outflowing gas). Filippenko et al. 
(1988) observe a remarkable double peaked structure for 
the Paschen lines, with speeds close to 300 km s -1 . 

4.3. Estimates of M and the windspeed 

If we assume that the observed radio emission is due 
to bremsstrahlung, we can obtain a rough estimate for the 
mass-loss rate, M, in this equatorial outflow. For this pur- 
pose, we assume that the outflow is radial but confined to 
an angle a with respect to the orbital plane of the binary. 
For a simple wind mass-loss law, the mass density, p, of 
the outflow then depends on the distance r from the sys- 
tem according to p = M / (47r sin a r 2 Voo), where Woo is the 
outflow velocity at infinity. At a particular radio frequency 
v, the outflow will be optically thick to some distance f„. 

Assuming that we see all the radio emission from the 
optically thin part of the outflow and are observing the 
system close to the orbital plane (both assumptions are 
only approximately true and ignore geometrical complica- 
tions) , a rough estimate for the mass-outflow rate is 

M ~ 1.6 x 10~ 4 Mq yr -1 (1) 

„3/4 ,3/2 -1/2 --1/2 / . xl/4 

xS 50 d 3 ^300^1.4 9w (sina) 3 (, , 

where 6*50 = 5^/50 mJy, d 3 = d/3kpc, u 30 o = 
Uoo/300kms _1 , V1.4 = j//1.4GHz, 510 = 5/IO (g is the 
Gaunt factor for free- free emission; see e.g., Rybicki & 
Lightman 1979), (sin 01)30 — sin a/ sin 30°. 

One of the major uncertainties in this estimate is the 
velocity of the outflow, though a velocity of ~ 300 km s -1 
is similar to that of the lines by Filippenko et al. (1988) 
and is close to the characteristic orbital velocity of SS 433, 
as one might expect for an outflow from the binary sys- 
tem rather than either binary component. Furthermore, 
if this outflow started soon after the supernova explosion 
which formed the compact object ~ 10 4 yr ago and whose 
impressively circular remnant is seen clearly in the images 
of Dubner et al. (1988), a velocity of ~ 300 kins" 1 implies 
an extent of the outflow of ~ 3arcmin (for d = 3kpc). 
Indeed, this is exactly the size of the extended smooth 
emission seen by Dubner et al. (1998) and suggests that 
this may be the outer extent of the same outflow. 

The above estimate for M would imply that the outflow 
is optically thick at a frequency of ~ 1 GHz to a distance of 
r ~ 10 15 cm and that the radio emission from the central 
region of the jet would be somewhat attenuated (although 
this will also depend on the exact geometry of the outflow 
and its orientation to the line of sight). 

It also suggests that the outflow is moderately optically 
thick in the optical and that part of the observed visual 
extinction to the system (Ay = 7.8, Margon 1984) may be 
due to the outflow, as postulated by Zwitter et al. (1991). 

The inferred mass- loss rate, M ~ 1O~ 4 M yr 1 , is 
much higher than any reasonable mass-loss rate from an 
O-star primary and suggests that it is connected with the 
unusual short-lived phase SS 433 is experiencing (see § 4.1). 
It could be mass loss from a common envelope that has al- 
ready started to form around the binary, or a hot coronal 
wind from the outer parts of the accretion disk driven, e.g., 
by the X-ray irradiation from the central compact source. 
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5. SUMMARY 

With unique sampling in the UV-plane, we have im- 
aged the SS433 system at 6 cm and 20 cm and securely 
detected at both wavelengths smooth emission extending 
over a few hundred AU perpendicular to the jet axis. The 
most likely interpretation of this radiation is emission from 
matter which has been ejected from the disk as a thermal 



wind with an outward speed of ~ 300 km s . 
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Frequency [GHz] Perpendicular distance from jet-axis [mas] 

Fig. 3. — a: The total flux density in the 'rufF emission as a function of frequency (see text). Crosses: northern emission; triangles: 
southern emission; filled circles: sum of northern and southern emission, b: The flux density integrated over 40 mas strips parallel to the jet, 
as a function of distance perpendicular to the jet, at 10 mas resolution. The solid line is 18 cm, the dotted line the 6 cm data. Note the flat 
spectrum of the 'ruff' emission, compared to the inverted spectrum of the (self-absorbed) jet core. 



